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Mass and isospin dependence of symmetry energy coefficients a sym of finite nuclei is investigated 
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with the measured nuclear masses incorporating the liquid drop mass formula. The enhanced 
o sym for nearly symmetric nuclei are observed. To describe the mass and isospin dependence of 
a symj a modified formula based on the conventional surface-symmetry term is proposed and the 
corresponding rms deviation of nuclear masses is checked. 
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The study of nuclear symmetry energy has attracted great attention in recent years both 
theoretically and experimentally llr]5|. The mass dependence of symmetry energy coefficients 
has been clearly observed and the obtained symmetry energy coefficients of finite nuclei are 
considerably smaller than the symmetry energy coefficient of nuclear matter at saturation 
density. The symmetry energy coefficient of finite nuclei is usually extracted by directly 
fitting the measured nuclear masses with different versions of the liquid drop mass formula. 
Some different forms for describing the mass dependence of symmetry energy coefficients of 
finite nuclei, which divide the symmetry energy of a nucleus into the volume and surface 



contributions, have been proposed in Refs. 



-|5[. The volume symmetry energy of nucleus 



corresponds to that of nuclear matter at saturation density. In this work, we investigate the 
symmetry energy coefficients of nuclei, especially the isospin dependence of the symmetry 
energy coefficient, based on the more than 2000 precisely measured nuclear masses J6(. 

We start with the well known liquid drop formula. The liquid drop energy of a nucleus 
is described by a Bethe-Weizsacker mass formula |7J], 

E LD (A, Z) = a v A + a s A 2/3 + E Coul + a sym I 2 A } (1) 

neglecting the pairing term. Where, I = (N — Z)/A denote the isospin-asymmetry. The 
symmetry energy coefficient a sym is conventionally expressed as a function of mass number 
A [3j. The Coulomb energy is written as 

£ Cou i = *c Z(Z A ~ l) (1 - 0.76Z- 2 / 3 ) (2) 

with the coefficient a c = 0.71 MeV. Inserting Eq.(2) into Eq.(l) and using the relation 
Z = ^(1 — /), the liquid drop energy per particle £ld = Eix>/A can be expressed as a 
function of mass number A and isospin asymmetry /. Performing a partial derivative of 
£le>(A -0 with respect to the isospin asymmetry /, the symmetry energy coefficient can be 
expressed as 
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ro) _ f de LD {A,I) l dE Coul {A,I) \ 
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FIG. 1: (Color online) Symmetry-energy coefficients of nuclei as a function of nuclear mass number 
(a) and of isospin asymmetry (b) . The shades denote the results of Danielewicz et al [3( . The crosses 
and the short dashes denote the extracted as ym and ag ym terms of the symmetry energy coefficients 
from the measured nuclear masses, respectively. 



Omitting the microscopic shell and pairing corrections and the corrections from nuclear 
deformation, the values of £LO dI ' could be obtained from the measured energy per particle 



-exp 



6], 



de LB (A, I) e^JA, I 2 ) - e cxp (A, I 



(5) 



dl I 2 - h 

Where, I\ and I-i denote the isospin asymmetry of nuclei (A,Z — AZ) and (A,Z + AZ), 
respectively. We take AZ = 1 in this work. Through iterations to a sym in Eq.(3), 
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one can obtain the expression of asym which includes all terms with 
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in Eq.(6), 
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In Fig.l (a), we show the extracted symmetry energy coefficients of nuclei as a function of 
nuclear mass number. The crosses and the short dashes denote the extracted as ym and ai ym 
terms of the symmetry energy coefficients from the measured nuclear masses, respectively. 



One can see that the contribution of Osym term is much smaller than that of ai ym term for 
most nuclei since I is a small quantity. The relatively large fluctuations in ai ym for heavy 
nuclei are mainly caused by the shell effects. The shades denote the results of Danielewicz 
et al. 
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The extracted ai ym term of the symmetry energy coefficients for heavy nuclei are comparable 
to the results of Danielewicz et al. For intermediate and light nuclei, there exists obvious 
oscillations and fluctuations in the extracted Osym which are probably caused by the shell 
effects and other nuclear structure effects. In the region A < 120, the extracted a sym are 
generally higher than the results of Danielewicz et al. In Fig. 1(b), we show the same data as 
in Fig. 1(a), but as a function of isospin asymmetry /. One can find that the obtained a^ym 
term of the symmetry energy coefficients somewhat depend on the corresponding isospin 
asymmetry of nuclei, especially for nearly symmetric nuclei. The a^ym term of the symmetry 
energy coefficient obviously increases with the decrease of asymmetry. The dependence of 
symmetry energy coefficient on the asymmetry of nucleus, especially that a sym increases with 
increasing proton fraction of the system is also found in 1|. The results from Danielewicz 
et al. cannot reproduce the observed trend of isospin dependence well. 

To describe the isospin and mass dependence of symmetry energy coefficient of nucleus, 
we propose a modified formula, 

2- \I\ 
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(9) 



A 1 / 3 ' 2+|/|A_ 

based on the conventional surface-symmetry term of liquid drop model, with a small correc- 
tion term from isospin asymmetry. The introduced correction term approximately describes 
the Wigner effect |8l-ll2| of nuclei. The introduced / term in a sym roughly leads to a correction 
Ew to the binding energy of the nucleus, 

" 2-\I\ ' 

_2+\I\A_ 

which is known as the Wigner term. In Fig. 2, we show the Wigner energies of nuclei 
calculated with different models. The open circles and the crosses denotes the results in 
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Ref . [lOj and those of this work, respectively. The straight line denotes the results of Satula 



et al. 11| . i.e. E\y w 47|J|. In 12], Myers and Swiatecki write the Wigner term as 
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FIG. 2: (Color online) Wigner energies of nuclei calculated with different mode 



and the straight line denote the results in Ref. 



10] 



s. The open circles 



111 ], respectively. The 



and those of Satula et al. 

crosses denote results of this work, with c sym = 29.3 MeV determined by fitting the 2149 measured 
nuclear masses. 

E w = -C exp[-W\I\/C } « -C + W\I\ + ..., with Cn = 10 MeV, W = 42 MeV. The 



results of this work are comparable to those from Ref. 11] for most nuclei. 

With increasing of mass number A, the a sym in Eq.(9) has a finite value approaches c sym 
which corresponds to the symmetry energy coefficient of nuclear matter at saturation density. 
The results from Eq.(9) are shown in Fig.3 for comparison, with c sym = 31 MeV and k = 2 



,(°) 



(o) 



which are obtained by fitting the extracted as ym . One can see that the extracted Osym can be 
reproduced reasonably well. Furthermore, we have checked the rms deviations of 2149 masses 
of nuclei with iV and Z > 8 from the measured data defined as a 2 = — ^ ( Aftxp — M^ ) 
by taking different forms of symmetry energy coefficients mentioned above incorporating 
the liquid drop mass formula of Eq.(l). The obtained rms deviations and the corresponding 
parameters of the liquid drop formula are listed in Table 1. Adopting the form in Eq.(8), we 
obtain an rms deviation of 2.71 MeV. With Eq.(9) for the symmetry energy coefficient, the 
rms deviation is reduced to 2.55 MeV. Compared with the case without the I term being 
taken into account, the rms deviation is reduced by 6% (see Table 1). Incorporating the 
semi-empirical mass formula in [13j], the rms deviation of the 2149 masses of nuclei can be 
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FIG. 3: (Color online) The same as Fig.l, but with the results of Eq.(9) (open circles) for compar- 



ison. 



TABLE I: rms a deviations between 2149 measured data and predictions of Eq.(l) with different 
a sym forms, and the corresponding optimal parameters of liquid drop formula. 

a sym form a v (MeV) a s (MeV) a c (MeV) c sym (MeV) re a (MeV) 



Danielewicz 


-15.55 


18.18 


0.71 


27.39 


1.28 


2.71 


c sym[l — reA J 


-15.57 


18.25 


0.71 


26.09 


0.80 


2.72 


This work 


-15.56 


18.11 


0.71 


29.38 


1.52 


2.55 



considerably reduced, falling to 0.516 MeV. Furthermore, when the isospin dependence of 
symmetry energy coefficient is taken into account, the obtained optimal c sym changes from 
26.09 to 29.38 MeV which is close to the calculated symmetry energy coefficient of nuclear 
matter at saturation density from the Skyrme energy density functional 13]. 



In summary, the mass and isospin dependence of symmetry energy coefficients a sym of 
finite nuclei has been investigated with the measured nuclear masses incorporating the liquid 
drop formula. For heavy nuclei, the extracted Osym term of the symmetry energy coefficients 
are consistent with the results of Danielewicz et al. For light and intermediate nuclei, 
there exists oscillations and fluctuations in the extracted asym- The isospin dependence 
of symmetry energy coefficients, especially the enhanced a sym in nearly symmetric nuclei, 
has been observed. To describe the mass and isospin dependence of a sym , we propose a 



modified formula based on the conventional surface-symmetry term, with which the isospin 
dependence of a sym can be described reasonably well and the rms deviation of nuclear masses 
from the experimental data can be effectively reduced. 

ACKNOWLEDGEMENTS 

This work was supported by National Natural Science Foundation of China, Nos 10875031 
and 10847004. 



[1] S. K. Samaddar, J. N. De, et al., Phys. Rev. C 76 (2007) 041602. 

[2] V. M. Kolomietz and A. I. Sanzlmr, Phys. Rev. C 81, 024324 (2010) 

[3] P. Danielewicz, et al., Nucl. Phys. A 818 (2009) 36. 

[4] M. W. Kirson, Nucl. Phys. A 798 (2008) 29. 

[5] J. Mendoza-Temis, J. G.Hirsch, A. P. Zuker, arXiv: nucl-th/0912.0882. 

[6] G. Audi, A.H. Wapstra and C. Thibault, Nucl. Phys. A 729 (2003) 337. 

[7] H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 82 (1936). 

[8] P. Moller, J. R. Nix, et al., At. Data and Nucl. Data Tables 59 (1995) 185. 

[9] D. Lunney, J.M. Pearson, C. Thibault, Rev. Mod. Phys. 75 (2003) 1021. 

[10] G. Royer, Roman. Rep. Phys. 59 (2007) 625. 

[11] W. Satula , D. J. Dean, et al., Phys. Lett. B 407 (1997) 103. 

[12] W.D. Myers, W.J. Swiatecki, Nucl. Phys. A 612 (1997) 249. 

[13] Ning Wang, Min Liu and Xizhen Wu, Phys. Rev. C 81 (2010) 044322. 



